Inconsistencies of the Adiabatic Theorem and the Berry Phase 
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The adiabatic theorem states that if we prepare a quantum system in one of the instantaneous 
eigenstates then the quantum number is an adiabatic invariant and the state at a later time is 
equivalent to the instantaneous eigenstate at that time apart from phase factors. Recently, Marzlin 
and Sanders have pointed out that this could lead to apparent violation of unitarity. We resolve the 
Marzlin-Sanders inconsistency within the quantum adiabatic theorem. Yet, our resolution points to 
another inconsistency, namely, that the cyclic as well as non-cyclic adiabatic Berry phases may vanish 
under strict adiabatic condition. We resolve this inconsistency and develop an unitary operator 
decomposition method to argue for the validity of the adiabatic approximation. 
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I. INTRODUCTION 

Adiabatic theorem is one of the most important and 
widely studied theorems in quantum mechanics [1,2]. It 
states that if we have a slowly changing Hamiltonian that 
depends on time, and the system is prepared in one of the 
instantaneous eigenstates of the Hamiltonian then the 
state of the system at any time is given by an instanta- 
neous eigenfunction of the Hamiltonian up to multiplica- 
tive phase factors [3-5]. It has potential application in 
diverse areas of physics such as in molecular physics, nu- 
clear physics, chemical physics, quantum field theory and 
so on. It is the revisit of adiabatic theorem that enabled 
Berry to discover his now famous geometric phase [6] 
which has numerous applications ranging from quantum 
Hall-effect [7] to understanding of damping of collective 
excitations in finite fermi systems [8]. Furthermore, adi- 
abatic theorem and Berry phase have played important 
roles in design of quantum algorithms [9] and implemen- 
tation of geometric quantum computations [10]. 

Recently, Marzlin and Sanders (MS) have pointed out 
an inconsistency in the quantum adiabatic theorem [11]. 
In this paper we resolve their inconsistency. However, our 
resolution leads to another inconsistency, namely, under 
strict adiabatic approximation cyclic as well as non-cyclic 
Berry phase almost vanish! We suggest that while deal- 
ing with adiabatically changing Hamiltonians one should 
not be within the strict regime nor one should be too 
much away from adiabatic regime. One has to optimize 
the operational scale so as not to encounter MS type in- 
consistency or vanishing adiabatic Berry phase type in- 
consistency. The organization of our paper is as follows. 
In section II, we discuss the Marzlin-Sanders inconsis- 
tency. In section III, we give a simple resolution to it. 
In section IV, we point out an apparant inconsistency 
with the cyclic and non-cyclic Berry phase and provide 
a resolution to it. In section V, we develop a unitary op- 
erator decomposition which can be split into a diagonal 



and off-diagonal part that captures the adiabatic theo- 
rem, Berry phase and transition to non-adiabatic regime 
in an elegant way. Also, we derive a lower bound for 
the smallness parameter that plays a crucial role in va- 
lidity of the adiabatic theorem. In section VI, we provide 
two examples to illustrate our novel results. Finally our 
conclusion follows. 



II. THE MARZLIN-SANDERS INCONSISTENCY 

Consider a quantum system that is governed by a 
Hamiltonian H(t) with discrete, non-degenerate spec- 
trum {E n (t)} (n — 0, 1, 2, . . .N). The instantaneous 
eigenstates \n(t)) of the Hamiltonian satisfies an eigen- 
value equation H(t)\n(t)) = E„(t)\n(t)) . If the state of 
the system was initially in one of the eigenstate (say nth), 
then under adiabatic approximation (AA) the state at a 
later time t is given by 



IV>(*)> 



*(*) e *7»W| n ( t )) 



(1) 



where 8 n (t) = — J E„(t')dt' (h = 1) is the dynamical 
phase and 7„(i) = / i(n(t')\h(t'))dt' is the extra phase 
that is usually neglected in the text book but can give 
rise to Berry phase during a cyclic change of the Hamilto- 
nian [6]. (Some authors mention that "f n (t) is the Berry 
phase which is incorrect. Unless we have a cyclic evolu- 
tion 7„(£) is not gauage invariant, hence not observable.) 

Since the MS inconsistency is a very curious one, (as 
well as for the sake of completeness), we spell out their 
argument. Let us define a state |V>(£)) = U ' (t)\t/j(0)) , 
where U(t) = Texp(— i J H(t) dt). This state satisfies a 
Schrodinger equation with a Hamiltonian H(t), i.e., 



(2) 



where H = -W(t)H(t)U(t). Note that This holds irre- 
spective of adiabatic approximation. Now, what is the 
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solution to the above equation in the adiabatic approxi- 
mation? First, note that if \n(t)) is an eigenstate of H(t) 
with eigenvalue E n (t), then \n(i)) = W(t)\n(t)) is an 
eigenstate of H(t) with eigenvalue —E n (t). Since |n(t))'s 
are orthonormal and satisfies completeness criterion, this 
forms an instantaneous basis. We can expand 



|#)} = ^d„(i)K*)), 



(3) 



and using A A the amplitudes d n (t)'s are given by 

id n {t) + E n {t)d n {t) + i(n(t)\n(t))d n (t) 

+ ^i{n(t)\ih(t))d rn (t)=0. (4) 



Under adiabatic approximation we can drop the terms 
like (fi(t)\m(t)) because 



\W)\Mt)) 

\E n — E m \ 



« 1, n =/= m. 



(5) 



In general if H(t) varies slowly a unitarily related Hamil- 
tonian need not vary slowly. But in this case it happens 
to be so. Under the above condition the solution is found 
to be 

d n (i)«exp(-i5 n (i))exp(i J i(n(t)\n(t)) dt)d n (0). (6) 

Therefore, the state |V>(t)) is given by 

$(t)) wexp(-i<y n (t))exp(i J i(n(t)\n(t)) dt)\n{t)) . (7) 

Using \n(t)) = exp(-iS n (t)) exp(-i^ n (t))\n(0)) and 
i(n(t)\n(t)) — i{n{t)\h(t)) — E n {t) one will have 



J E n (t')dt' 



n(0)). 



(8) 



This was the result of MS. Now we give their contradic- 
tion. Since |V>(t)) is unitarily related to the initial state, 
so it must be normalized to unity. However, the above 
solution together with the standard adiabatic ansatz (1) 
gives a non-unit norm! Explicitly, one can see that 

(m\m) = (m\u(t)\m) « e ^«<n(o)K*)) ^ 1 

(9) 

Hence a contradiction. So the standard adiabatic theo- 
rem apparently violates unitarity! Does it really do so? 



III. RESOLUTION OF THE MS 
INCONSISTENCY 

The resolution to the above contradiction is now given 
within the adiabatic theorem. Now let us take a close 



look at the transition amplitude A n (t) = (n(0)\n(t)) be- 
tween the initial and the instantaneous eigenstates, and 
ask how does that change with time as we slowly change 
the Hamiltonian. Consider the transition amplitude de- 
fined by 

A n (t) = (n(0)|n(t)>. (10) 
Its time rate of change is given by A n (t) 
Mt) = =5>(0)|m(t))<m(t)|n(t)) = (n(0)\n(t))(n(t)\n(t)) 

m 

+ ^(n(0)\m(t))(m(t)\h(t)). (11) 

Now under standard adiabatic approximation one can 
drop the terms (m(t)\n(t)), for m ^ n. Then we have 
the transition amplitude A n (t) as follows: 



iA n {t) w j n {t)A n {t), 



(12) 



where j n (t) — i(n(t)\h(t)) is the Berry frequency. This 
leads to 



A n (t) w e-^A n (0). 



(13) 



Since A n (0) = 1, this implies that (n(0)\n(t)) w e" l7 "W. 
Using this solution one can easily see that the MS contra- 
diction is resolved, i.e., the unitarity is preserved. This 
also tells us that the transition probability between initial 
eigenstate and the later instantaneous eigenstate is unity 
for all time. In terms of the 'minimum-normed distance' 
[12] we have L» 2 (|n(0)), |n(t))) - 2(1 - \(n(0)\n(t)}\) w 
which is almost zero. So under strict adiabatic condi- 
tion the instantaneous eigenstate apparently stays almost 
close to the original one and hence there is no violation 
of norm preservation. 



IV. THE BERRY PHASE- YET ANOTHER 
INCONSISTENCY 

The above resolution could have been a satisfying sit- 
uation if the following is not true. However, as we will 
show below under strict adiabatic evolution, i.e., under 
the condition (5) the cyclic as well as the non-cyclic Berry 
phases almost vanish! 

Consider the cyclic variation of the Hamiltonian, i.e., 
H(R(T)) = H(R(0)) over a period of time T. Then we 
know that the state of the system at time t — T is given 

by 

|V>(T)) w e^Ia En{t)dt e^ c )\iP(0)) (14) 
where 7 n (C) is given by 

7n (C)=»jf (n(t)\h(t)) dt = ^<n(R)|Vn(R)).dR (15) 
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is the gauge-invariant Berry phase that depends only 
on the geometry of the path in the parameter 
space and is also measurable [6]. Therefore, we 
have (V(0)|V(T)} = exp[i(5 n (T) + y n (C))]. How- 
ever, our solution to MS inconsistency suggest that 
(V(0)|VCO> = exp^„(T))exp( l7 „(C))(n(0)KT)) = 
exp(iS n (T)). This implies that the initial and the final 
state differ only by the dynamical phase and there is no 
observable Berry phase. Hence, a contradiction! 

Next, consider the general definition of the Berry phase 
when a pure state vector undergoes a time evolution 
| -0(0)) — ► \4>(t)). Invoking Pancharatnam's idea of rel- 
ative phase shift one can define the geometric phase for 
non-cyclic evolution of quantum systems [13]. The non- 
cyclic geometric phase is then given by a gauge invariant 
functional of ip(t) along an open path T [14] 



$ G [T] = ArgMOM*)) + 



i J dtm)\m) 



(16) 



Using the reference-section state vector \x(t)) we can ex- 
press the general geometric phase during an arbitrary 
evolution as 

<MT]=t J dt( X (t)\x(t)) = I G, (17) 

where the reference-section \x(t)) — |(^[t)|^(o))| IV'ffl) an< ^ 
G = i(x\dx) is the generalized gauge potential or connec- 
tion form [15,16]. Thus, the generalized geometric phase 
can be written as a line integral of a vector potential in 
the projective Hilbert space of a quantum system. For 
differential geometric formulation of the general Berry 
phase see [15]. 

Note that the above definition holds irrespective of adi- 
abatic, cyclic, and Schrodinger time evolution. So this 
is the generalized geometric phase during a time evo- 
lution of a quantum system described by a pure state 
vector. Under adiabatic approximation there is an open- 
path Berry phase which was introduced in [17] and is 
given by 



[r] = Arg<n(0)|n(i)) + i J dt(n(t)\h(t)) 
= ^G„(R).dR, 



(18) 



where G n (R) is the generalized gauge potential that 
gives rise to the adiabatic open-path Berry phase. Under 
strict adiabatic approximation using our previous solu- 
tion we find that [r] « 0. Thus, it almost vanishes for 
all time! Therefore, it appears that under strict adiabatic 
evolution a quantum system cannot acquire any Berry 
phase (non-cyclic as well as cyclic). However, there are 
many physical systems that show the existence of non- 
cyclic and cyclic Berry phases under adiabatic approx- 
imation. Thus our solution , though resolves Marzlin- 
Sandcrs inconsistency, yet points to another important 
inconsistency. 



Hence, a possible way out is not to drop off-diagonal 
terms that are usually done. One has to be careful when 
to drop and when not to. The statement that matrix ele- 
ments causing transition to other cigenstates are dropped 
under adiabatic approximation could lead to internal in- 
consistencies either MS type or our type (i.e. the vanish- 
ing Berry phase) . Thus the source of our inconsistency 
is dropping of the terms like (m(t)\h(t)) , for m ^ n. 

We have seen that under strict adiabatic approxima- 
tion the transition amplitude between the initial and 
the instantaneous eigenstate obeys a linear homogeneous 
equation that results in a vanishing Berry phase. How- 
ever, if we investigate the full solution and one can save 
the adiabatic Berry phase. Note that Eq(4) can be writ- 
ten as 



ifcAni 1 ) =in{t)A n {t) 



Sn(t), 



(19) 



where S n (t) = £^<n(0)|ro(t))<m(t)|n(i)). The solu- 
tion to the above equation can be written as 



A n (t) 



-iln(i) 



[1 



dt'S n (t')e 



(20) 



The second term clearly represents the correction to the 
adiabatic approximation and its presence can only make 
a non-zero Berry phase in cyclic as well as non-cyclic 
case. In particular, the generalized non-cyclic adiabatic 
Berry phase given by Eq(8) can be expressed as 



4">[r] 



tan 



■'ti 



-ReQ n {t) , 
+ ImQ n {ty 



(21) 



where Q n (t) = J dt' S n (t')e l ^ (t '\ This can be shown 
to be related to the response function of a many body 
quantum system that explains damping of collective ex- 
citations [8]. 



V. UNITARY OPERATOR FOR ADIABATIC 
EVOLUTION 



We can develop a general solution to the unitary evolu- 
tion operator and show that it has two pieces; a diagonal 
and a non-diagonal piece. It is the diagonal piece that 
gives what we want - adiabatic theorem but there is no 
inconsistency if we keep the order of approximation in 
our development. 

The time evolution operator of a quantum system 
obeys 



iU(t) = H{t)U{t) 



(22) 



with U (0) = /, and the unitarity condition for all t holds, 
i.e., U{t)U{ty = U(tfU(t) = I, where / is the idenity 
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operator. Since the Hamiltonian is Hermitian for any 
t, it admits an eigen-expansion (for simplicity, we con- 
sider here discrete and non-degenerate case) H(t)\n(t)) = 
E n {t)\n(t)), (n(t)\m(t)) = 5 nm and £ n \n{t)){n{t)\ = 7, 
for all t. We can express the unitary operator in terms 
of these instantaneous eigenstates as 



U{t) = Y,Ht))U nm {t){m{% 



(23) 



where the matrix elements satisfy 

iU nm (t) + E n (t)U nm {t) + 5^i<n(t)|p(t))^ pm (t) - 0. (24) 



The condition at t = implies that U nm (0) = 
Snm = ^mn(O) an d the unitarity relation implies that 
S P ^pn(i)^m P (*) — <5«ni f° r au Now let us introduce a 
"smallness" parameter, e, such that 



U nm {t) = U nn {t)5 nrn + e5U nm {t){\ - S nm ). 



(25) 



Then the unitarity condition to leading order in e, im- 
plies \U nn (t)\ 2 = 1 for all t. Hence, U nn {t) = e^"W, 
with [/„„(0) = 1 or </> n (0) = 0. The explicit form of the 
phase to this order is given by the real number 



.(*) 



t rt 

E n {t')dt' + i(n{t')\h{t')) dt' . (26) 



The first term is the dynamical phase and the second one 
gives us the familiar Berry term. The corresponding U(t) 
is then of the form 



U(t)=J2^ {t) Ht))(n(0)\ 



+ eJ2Ht))SU nm (t)(m(0\. (27) 

Thus the adiabatic theorem is presented with a consistent 
form. 

The 'smallness' parameter actually decides when to 
keep the off-diagonal terms and when the approxima- 
tion with diagonal ones is satisfactory. But how small 
can it be? One can obtain a non-trivial lower bound 
on the e as follows. Let the initial state of the sys- 
tem is IV'(O)) = \n(0)). The state at a later time t is 
\tp(t)) = U(t)\n(0)). Now the transition amplitude be- 
tween the initial and the final state is 

(m\m) = (n(0)\U(t)\n(0)) = U nn (t)(n(0\n(t)) 

+ €^2SU mn {t)(n{0\m{t)). (28) 



This leads to 



\(m\m)\ < \(n(0\n(t))\+eJ2 (29) 



The above inequality can be written as 

D(\n(0),\n(t)))-D(\m),\m)) 



e > 



\SU mn (t)\ 



(30) 



where D(\n(0), \n(t))) is the 'minimum-normed' distance 
functions between the initial eigenstate and final eigen- 
state, and D(\ip(0)) ,\ip(t)}) is the similar one between 
the initial state and the final state of the system [12]. 
When we are within adiabatic regime D(\n(0), \n(t))) and 
D(\tp(0)}, \ip(t))) are same, hence we have e = 0. There is 
no lower bound on it. However, if we move away from adi- 
abatic regime D(\n(0),\n(t))) and D(\ip(0)),\ip(t))) dif- 
fer. Then we will have a lower bound on the 'smallness' 
parameter. 

One can also resolve the MS inconsistency us- 
ing the unitary operator method. Note that we 
can obtain \n(t)) = £ p b(0))t/ p *„ = U* n \n(0)) + 
off-diagonal terms, and thus as before we have 
\n(t)} w exp(i J E n (t')dt') exp(-ry„(i))|n(0)). Then us- 
ing (n(0)\n(t)} fa exp(— Z7„(i)) one can resolve the MS 
inconsistency. Below we give two examples to illustrate 
the power of unitary operator method developed in this 
letter. 



VI. TWO EXAMPLES 

Here we consider two examples that will illustrate the 
main points. 



A. The MS Example in the instantaneous 
representation 

Let us consider the precession of a spin-half particle 
in a magnetic field of strength proportional to u>o and in 
addition it rotates in the x — y plane with a frequency 
ft = 2tt/t. This is also considered by MS [11]. We work 
here in the instantaneous representation. The Hamilto- 
nian is given by H(t) = R(i).<r, where 



§(l-cos2w i) v~>, - 2 — 

1 ; ^ e int (w + i § sin 2w i) - § ( 1 - cos 2cj t) 

Its eigenvalues and eigenfunctions are E\(t) 

w 2 + ft 2 sin 2 = Eo(t), £?i(0) = w ; |m(t)) 
col(ai(t),6i(t)), |m(0)) = l/V2col(l,l), with 



e- tnt (u> -i§sin2w 



„(+\_ I ^oft) + ^ sin2 ^0* o -i(nt+9(t))/2 



h(t) = 



j E Q (t) - sin 2 u Q t ^ (nt+f)(t))/2 
2E (t) 



e™~. (31) 
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And similarly, we have E 2 (t) = —E (t), E 2 (0) = -ujq; 
\n 2 {t)) = col(a 2 (t),& 2 (t)), |n 2 (0)) = l/>/2col(-l, 1), 
with 



02 (*) 



'£<)(<) - ftsin 2 w t 



2£ (i) 



-»(nt+e(t))/2 



2E (t) 



(32) 



In the above equations, 9(t) = tan _1 [(fi/2a;o) sin2u;ot]. 
The unitary time evolution associated with the Hamilto- 
nian H(t) is given by 

cosuo t — ie~ tQt sinui t " 

COS LU()t 



U(t) = 



-ie lQ,t sin u t 



a y sm6sm(f>+a z cos6) = -gnoH { sindei<t> 



We need to express this in terms of the instantaneous 
eigenstates to isolate the adiabatic term from the non- 
adiabatic pieces. We present here the diagonal terms, 
being the adiabatic contributions to the evolutions: 
(m(t)\U(t)\ni(0)) = ^[ai(t)(cosw ^ - ie" int sin uj t) + 
bi(t)(cosuj a t - ie int smoj t) and (n 2 (t)\U(t)\n 2 (0)) = 
-j=[a 2 (t) (cos uJot + ie~ tnt sin uj n t) — b 2 (t) (cos iv t + 

ie sin cuot) ■ These exact solutions can be evaluated 
explicitly and compared with the solutions based on 
adiabatic approximation. 



B. The Schwinger example 

This example is an elegant one originally due to 
Schwinger [5] for describing the precession of a spin in 
a transverse time-dependent field. Here the Hamiltonian 
of the system is given by H(t) = —g^H(a x sin 9 cos <f) + 

cos 9 sin9e~^' 
- cos 9 

where the magnetic field H and 9 are independent 
of time, and (f> = u)t. The MS model is a version 
of the Schwinger's NMR precession problem described 
above. The instantaneous eigenvalues are E\ — g^oH 
and E 2 = —gfi H. The instantaneous eigenstates are 
|m(i)) = col(e-^/ 2 sin (9/2, -e"^ 2 cos 9/2) and \n 2 (t)) = 
col(e-^/ 2 cos(9/2,e^/ 2 sin(9/2). From these we have 
(m(t)\hi(t)) = i(uj cos9)/2 = -(n 2 (t)\h 2 (t)). Similarly, 
we have (ni(t)\n 2 (t)) = — i(u>sm6)/2 = (n 2 (t)\hi(t)) . 
The solution to the time evolution operator U(t) in 
terms of the instantaneous eigenstates can be obtained 
as U(t) = £ 2 . =12 MWij(t){nM\ ^th E^-(O) = 
The equations to be solved are The matrix elements of 
the time evolution operator obey Eq(13). 

iUu + i{ni\h 2 )U 2 i = (E 1 - i(ni\hi))Uu 

iU 21 + i(n 2 \hi)Un = (E 2 - i(n 2 \n 2 ))U 21 . (33) 

A similar pair of equations for the other two. We can 
solve these equations by Laplace transforms and obtain 
the solutions: The solutions are given by 



U u (t) = [Ei cosEit 

- i{gvoH + | cos9) smEit]/Ei = u; 2 (t) 

U 21 (t) = i[uj sin 9 sin E 1 t}/2Ei = U 12 (t), (34) 

where Ei = [(gfi H) 2 + g^ Q Huj cos9 + uj 2 /4} 1 / 2 . The 
unitarity of the this JJ-matrix is obeyed, as it should be. 
These are the exact solutions of the Schwinger problem. 

If we make the adiabatic approximation, we would 
drop the U 2 i term in Eq.(22) and obtain the following 
result: 



Un(t) w ex.p[-it(gfj, H + - cos 9)} = U£ 2 (t) 
U2l{t) " l 2 ((c/2) cos 9 + g^H) ] 



U 12 (t). (35) 



This means that we drop terms of the order of (to/2) sin 9 
in the exact expressions, so that E\ k, [g^oH + 
(uj/2) cos 9]. We thus see that the "smallness" parameter 
in the "adiabatic" treatment is (ni(t)\n 2 (t)). Note that 
the unitarity of this U (t) is verified to be obeyed con- 
sistent with the adiabatic approximation stated above. 
Also, a way to see the departure of the result in the adi- 
abatic approximation from the exact result is to com- 
pute the overlap A{ni(t)\ni(t)) e and hence the fidelity, 
F{t) = \A{ni{t)\ni(t)) e\ 2 ■ The overlap is given by 
A(ni(t)\ni(ij) E = exp(iflt)(cosEit- i(£l/ Ei) sin Eit) + 
(lu sin 9 /2) 2 (sin Erf /Ei) [sin nt/il), where = (g^H + 
(u)/2) cos 9). One sees at once that to the leading order 
in (u) sin 9/2), the fidelity is unity. 



VII. CONCLUSIONS 

Adiabatic theorem though widely studied, can reveal 
surprising phenomena, and sometime even apparent in- 
consistencies, namely, MS type or vanishing Berry phase 
type. Thus, it is of utmost importance to know where 
such inconsistencies may arise. Besides showing how to 
resolve these inconsistencies associated with the adia- 
batic theorem (AT), we have presented here several novel 
results. We have established an equation of motion for 
the transition amplitude leading to the Berry phase, de- 
veloped the unitary time evolution operator expressed in 
instantaneous basis as a powerful tool in AT as illustrated 
with two examples, estimated the "smallness" parame- 
ter in the unitary operator, and quantified the departure 
of AT from the exact by means of "Fidelity" expressed 
in terms of the respective unitary evolution operators. 
Given the current interest in the geometric phases for 
mixed states [18], this work opens up the possibility of 
studying the Berry phase and adiabatic theorem for den- 
sity operators and entangled quantum systems in a new 
perspective which one may like to take up in future. 
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Note Added: After completion of our work we noticed 
Ref. [19] which addresses the MS inconsistency. Our 
present work goes beyond this aspect of the adiabatic 
theorem. 
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